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We study the space-time evolution of "fiavor" neutrino wave-packets at finite temperature and 

density in the early Universe prior to BBN. We implement non-equilibrium field theory methods 

and linear response to study the space-time evolution directly from the effective Dirac equation in 

the medium. There is a rich hierarchy of time scales associated with transverse and longitudinal 

dispersion and coherence. A phenomenon of "freezing of coherence" is a result of a competition 

between longitudinal dispersion and the separation of wave-packets of propagating modes in the 

medium. Near a resonance the coherence and oscillation time scales are enhanced by a factor 

1/ sin 29 compared to the vacuum. CoUisional relaxation via charged and neutral currents occurs 

on time scales much shorter than the coherence time scale and for small vacuum mixing angle, 

shorter than the oscillation scale. Assuming that the momentum spread of the initial wave packet 

^^ . is determined by the large angle scattering mean free path of charged leptons, we find that the 

*vj ' transverse dispersion time scale is the shortest and is responsible for a large suppression in both the 

survival and transition probabilities on time scales much shorter than the Hubble time. For small 

C , mixing angle the oscillation time scale is longer than the collisional relaxation scale. The method 

^ ' also yields the evolution of right-handed wave packets. Corrections to the oscillation frequencies 

emerge from wave-packet structure as well as from the energy dependence of mixing angles in the 

C*^ , medium. 

(N 
> 

in . I- INTRODUCTION 

zi . Neutrinos are the bridge between astrophysics, cosmology, particle physics and nuclear physics [11,1^1^1^ ■ In recent 

f~^ ' years, there has been increasing experimental evidence that confirm that neutrinos are massive and oscillate between 
\Q .. different flavors[^la>l3>Bl3 providing the first indisputable hint of new physics beyond the Standard Model. 
f^ ■ Neutrino mixing and oscillations in extreme conditions of high temperature and density play a fundamental role in 

astrophysics and cosmology , Id Ull lia . [l3 . [l4 llSl [lq| . Resonant flavor mixing due to Mikheyev-Smirnov-Wolfenstein 
(MSW) effect can provide a concrete explanation to the solar neutrino problem [13, 13 ■ During Big Bang Nucle- 
I osynthesis (BBN) neutrino oscillations may result in corrections to the abundance of electron neutrinos[10|. This in 

^^. turn changes the neutron-to-proton ratio, affecting the mass fraction of "^He (see ref. _10J and references therein). 
Neutrino oscillations violate lepton number leading to the possibility that the cosmological baryon asymmetry may 
originate in the lepton sector [l9l l20l |21| . 

Neutrino propagation in a cold medium has been first studied by Wolfenstein|l7| who included the refractive 

index from electron neutrinos. Early studies focused on the neutrino dispersion relations and damping rates at the 

^ ' temperature limit relevant for stellar evolution or BBN [22 ■ This work has been extended to include charged leptons, 

Cd [ neutrinos and nucleons in the thermal medium [23 • The matter effects of neutrino oscillations in the early universe 

have been investigated in \wL 122 . |2J, [23 . More Recently, a non-equilibrium field theoretical description of neutrino 

oscillations in the early universe in the real-time formulation has been reported in [26j- 

Kayser|23 first pointed out subtle but important caveats in the vacuum oscillation formula obtained from the stan- 
dard plane wave treatment, which result from assuming a definite neutrino momentum for different mass eigenstates. 
He showed that knowledge of momentum allows experiments to distinguish different neutrino mass eigenstates, essen- 
tially destroying the oscillation pattern. He then proposed a wave-packet treatment of neutrino oscillations, in which 
the neutrino momentum is spread out. Since then, the wave-packet approach has been studied by many authors in 
both quantum mechanical 1, 2, 28, 29, 30, 31, 32] and field theoretical 33, 34, 35, 36, 37] frameworks, including the 
study of oscillations of neutrinos produced and detected in crystals [38J. 

The quantum mechanical approach usually refers to the intermediate wave-packet model in which each propagating 
mass eigenstate of neutrino is associated with a wave-packet [23 . This model eliminates some of the problems in the 
plane wave treatment although several conceptual questions remaining unsettled|28l]. See |37| and references therein 
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for detailed descriptions of these issues. A field theoretical approach is the external wave-packet model [33 in which 
the oscillating neutrino is represented by an internal line of a Feynman diagram, while the source and the detector are 
respectively described by in-coming and out-going wave-packets. A recent reviewJ23 presents the different approaches, 
summarizes their advantages and caveats and includes the dispersion of wave packets in the study. 

An important physical consequence of the wave-packet description of neutrino evolution is the concept of the 
coherence length beyond which neutrino oscillations vanish. A "flavor neutrino" wave-packet is a linear superposition 
of wave-packets of mass eigenstates. The different mass states entail a difference in the group velocity and an eventual 
separation of the wave-packets associated with mass eigenstates. This separation results in a progressive loss of 
coherence as overlaps between the wave packets dimishes. See for example[31| for an early explanation. In an actual 
source-detector experimental setup, the observation time is usually not measured and is commonly integrated out in 
a wave-packet treatment J31|. This leads to a localization term in the vacuum oscillation formula, which states that 
neutrino oscillations are suppressed if the spatial uncertainty is much larger than the oscillation length[3l|. 

The coherence of neutrino oscillations in matter has been studied within a geometrical representation in J30j , but 
the medium oscillation formula was not derived. While most of the studies focus on reproducing the standard vacuum 
oscillation formula, a consistent study of neutrino mixing and propagation in a medium in real time has not yet 
emerged. 

While in the vacuum the space-time propagation can be studied in the wave-packet approach by focusing on the 
space-time evolution of initially prepared single particle "flavor states" , the study of the space-time evolution in a 
medium at finite temperature and density requires a density matrix description. 

To the best of our knowledge, a full finite temperature field theoretical treatment of the space-time propagation of 
neutrino wave-packets in a medium including medium corrections and dispersion dynamics has not yet been offered. 
We consider this study an important aspect of the program to study the non-equilibrium evolution of neutrinos in 
the early Universe. Detailed studies have shown that neutrino oscillations and self-synchronization lead to flavor 
equilibration before BBN|33, l40l llU 142. l4a |. beginning at a temperature of T '--^ 30MeV with complete flavor 
equilibration among all flavors at T ~ 2MeV [4l|. If neutrinos are produced in the form of spatially localized 
wave-packets rather than extended plane waves before BBN, the two mass eigenstates separate progressively during 
propagation due to the small difference in group velocities. A significant amount of neutrino oscillations, which 
are crucial for "flavor equalization" requires that the two mass eigenstates overlap appreciably throughout their 
propagation, namely the coherence time scale should be sufficiently large to ensure that "flavor equilibration" through 
oscillations is effective. Therefore, it is important to pursue a full field theoretical study of neutrino wave-packet 
propagation in the medium directly in real time to determine the relevant time scales for coherence to be maintained 
and to identify the processes that contribute to its loss. 

A. The main goals of this article 

In this article we study the space-time evolution of neutrino wave packets in extreme environments at high temper- 
ature and density, conditions that prevail in the early Universe or during supernovae explosions. Our goals are the 
following: 

• i:) to provide a consistent and systematic non-equilibrium field-theoretical formulation to study the space-time 
evolution of initially prepared neutrino wave packets at finite temperature and density. This goal requires a 
treatment of the space-time evolution in terms of the density matrix, which goes beyond the usual treatment 
in terms of single particle states. To achieve this goal we implement a recently developed method[2g to study 
non-equilibrium aspects of neutrino propagation in a medium as an initial value problem in linear response. This 
method yields the effective Dirac equation of motion for the expectation value of the neutrino field induced by 
an external source. The effective Dirac equation in the medium includes self-energy contributions from charged 
and neutral currents up to one loop. 

• ii:) To assess the different time scales associated with wave-packet dispersion, coherence and oscillations including 
the medium effects, in particular near possible resonances in the in-medium mixing angles. This is achieved by 
solving the effective Dirac equation in the medium, which includes self-energy corrections, as an initial value 
problem. The initial wave-packet configuration is "prepared" by an external source term in linear response. 
This method also allows to assess corrections from the energy dependence of the mixing angles in the medium 
upon the wave-packet dynamics. 

• iii:) The space-time evolution of the initially prepared wave packet, including dispersive effects allows an assess- 
ment of the competition between the progressive loss of coherence in the wave-packet dynamics by the separation 



of mass eigenstates, collisional decoherence and cosmological expansion. While our study only includes the self- 
energy from charged and neutral currents up to one-loop, the final result allows us to include results available 
in the literature[l3 122, |2J, |23 to understand the effects of collisional decoherence and cosmological expansion 
when there is a separation of time scales. 

• iv): We focus our study within the context of early Universe cosmology, in particular in the temperature regime 
just prior to BBN where there is a possibility for resonant transitions [lOl 122 . |24L l25l l2q| . Of particular interest 
are the medium modifications of the dispersion relations, wave-packet dispersion, oscillation and coherence time 
scales in this temperature and energy regime. 

• v): A bonus of this field theoretical formulation is that it also allows to obtain the space-time evolution of right- 
handed neutrino wave packets. Although the amplitude for such wave packets is suppressed by M /k with M, k 
the typical mass and energy scale of neutrinos in the medium, and may not be relevant for neutrino processes 
in the early Universe, the method systematically yields this information. 

Since we study the propagation of neutrino wave-packets in a medium, aspects associated with the source-detector 
measurement processes are not well-defined or relevant in this case. Consequently, in contrast to most studies in the 
literature, we do not integrate in time as is the case for a description of experiments in the vacuum [3l|. Therefore, 
our study of propagation is both in space and time. 

B. Main results: 

Our main results are the following: 

• A systematic field theoretical formulation of the space-time dynamics of wave packets of massive and mixed 
neutrinos in terms of the effective Dirac equation including self-energy corrections. The space-time evolution is 
approached as an initial value problem via linear response with the full density matrix. 

• Wave packet evolution features characteristic time scales associated with transverse and longitudinal dispersion. 
The ratio of these scales is given by the enormous Lorentz dilation factor in the case of relativistic neutrinos. 
Neither of these scales receives substantial medium corrections. The shortest scale, associated with the transverse 
dispersion dominates the suppression of both the survival and transition probabilities. There is an interesting 
phenomenon of "coherence freezing" which results from the competition between longitudinal dispersion and 
the separation of the mass eigenstates. We find that there are medium as well as wave-packet modifications to 
the oscillation formula both for the oscillation frequency as well as the survival and transition probabilities. 

• There is a resonance for the mixing angle in the medium just prior to BBNJl^, ll^, 124. 1251 126| at an energy 
and temperature k ^T ^ 3.6 MeV for large mixing angle or fc ~ T ^ 7MeV for small mixing angle. Both the 
coherence and the oscillation time scales are enhanced by a factor l/sin26' near the resonance, where 9 is the 
vacuum mixing angle. This suggests a substantial increase both in the coherence and the oscillation scales for 
1 — 3 mixing, but not an appreciable modification for 1 — 2 mixing. 

• Assuming that the momentum spread of the initial wavepacket is determined by the large angle scattering mean 
free path of charged leptons[l|, we find that near the resonance the loss of coherence via charged and neutral 
current elastic scattering is faster than the loss of coherence by the separation of the mass eigenstates and occurs 
on a time scale much shorter than the Hubble time. However, the decoherence time scale is many orders of 
magnitude larger than the time scale for transverse dispersion which ultimately determines the suppression of 
the survival and transition probabilities. 

C. Outline: 

This article is organized as follows. In section m we obtain the effective Dirac equation of neutrino in a thermal 
medium implementing the methods of non-equilibrium field theory and linear response J2y|. In this section we obtain 
the in-medium dispersion relations and mixing angles. In section IIIII we develop the general formulation to study 
the space-time propagation of neutrino wave-packet. In this section we discuss the different time scales associated 
with dispersion, oscillations and coherence. In section Hvl we compare the different time scales with the Hubble and 
collisional relaxation time scales and discuss the impact of the different scales upon the space-time evolution of the 
neutrino wave-packets, coherence and oscillations. We present our conclusions in section Ivl 



II. EFFECTIVE DIRAC EQUATION OF NEUTRINOS IN A MEDIUM AND LINEAR RESPONSE 

The study of the evolution of neutrino wave packets in the vacuum typicaUy involves a description of the experimen- 
tal production and detection of these wave packets. We study the space-time evolution of wave packets in a medium 
as an initial value problem. This is achieved in linear response by coupling an external source term which induces an 
expectation value of the neutrino field in the density matrix, after this source is switched off the expectation value 
evolves in time. The propagation of this initial state is described by the effective Dirac equation in the medium, 
which includes the self-energy corrections. This is the familiar linear response approach to studying the evolution 
out of equilibrium in condensed matter systems. The correct framework to implement this program is the real-time 
formulation of field theory in terms of the closed-time-path integral |4J, lisL liq . 

We restrict our study to the case of two Dirac flavor neutrinos, while the formulation is general and can treat 1 — 2 
or 1 — 3 mixing on equal footing, for convenience we will refer to electron and muon neutrino mixing. Neutrino mixing 
and oscillations is implemented by adding to the Standard Model a Dirac mass matrix Mab which is off-diagonal in 
the flavor basis. For our discussion, the relevant part of the Lagrangian is given by 

C = Cl+Cl,+C%+CcC + CnC + Va^a + ^aVa, (H.l) 

where L^ is the free field neutrino Lagrangian density 

Ll^Va{l^5al>-Mab)vb {11.2) 

where a, 6 are flavor indices. L^/z ^-^^ the free field vector boson Lagrangian densities in the unitary gauge, namely 
C°w = -\ {d,,W+ - duW+) {d^'W- '^ - d'^W- '') + M^ W+ W- ^ , (IL3) 
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i {d^Z, - d,Z^) {d^Z-^ - d'^Z'^) + i. 



^i = -7 id^^Z, - d.Z^) (S^Z" - d^'Z^') + -Ml Z^ Z^ , (IL4) 



and the charged and neutral current interaction Lagrangian densities are given by 

LCC = ^ ^a l" L la W+ + Ia7^ L Va W"] , (IL5) 

v2 
'^^^ = ^TtLt ^- ^' ^ "" ^^ + 7« ^'^ (^r - 9t 7') /a Z^\ , (IL6) 

Zi COS yJ ijj 

where B^, is the Weinberg angle, L = (1 — 7^)/2 is the left-handed chiral projection operator, and (7^'^ are the vector 
and axial vector couplings for quarks and leptons. The label I stands for leptons and / generically for the fermion 
species with neutral current interactions. The external sources r]a,fja which couple to the neutrino fields depend 
explicitly on space and time and induce an expectation value whose time evolution is studied in linear response. 
For two Dirac flavor neutrinos, the mass matrix Mab is parametrized by 



M = "" "^ . (IL7) 

For the vacuum case, the elements rueei f^^i^i and rriep are related to the vacuum mixing angle 9 and masses of the 
propagating mass eigenstates Mi and M2 as follows 

ruee = cos^ e Ml + sin^ 6 M2 ; m^^ = sin^ 9 Mi + cos^ 9 M2 ; rue^ = -{Mi - M2) sin 9 cos 9 . (IL8) 

For later convenience and to establish contact with observable parameters we introduce the following quantities 

j^^Mi + M2 . SM'^Mi-Ml (IL9) 



The masses Mi,2 can be conveniently written in terms of these quantities as 



Ma - M 



, .n 1 SAP 

AM 



■ a ==1,2. (11.10) 



The current value for M obtained by WMAP |43| and the oscillation parameters from the combined fitting of the 
solar and KamLAND data are |4l| respectively: 

M w 0.25 {eV) ; {SM^^l ~ 7.9 x 10"^ (eV)'^ ; tan^ 6*12 w 0.40 . (11.11) 

For atmospheric neutrinos, analysis from SupcrKamiokande, CHOOZ and atmospheric neutrino data yield, 

\5M^^\ « (1.3 - 3.0) X 10^3 (eV)^ ; sin^ 6*13 < 0.067 (3cr) . (11.12) 

2 

This implies that |(5M^|/ M ^ 1, an almost degenerate hierarchy of neutrino masses. 

A. Linear response: 

The medium is described by an ensemble of states, and the description is in terms of a density matrix. Therefore 
the question of space-time evolution is more subtle, while in the vacuum one can consider preparing an initial single 
particle state and evolving it in time, such a consideration is not available in a medium, and the question of time 
evolution must be formulated differently, namely in terms of expectation values of the relevant operators in the density 
matrix. 

In equilibrium the neutrino field cannot have an expectation value in the density matrix. The usual method in many 
body theory to study the non-equilibrium evolution of single quasiparticle states is the method of linear response: 
an external source is coupled to the field which develops an expectation value in the density matrix induced by the 
source. The expectation value of this field obeys the equation of motion with medium corrections. Upon switching-off 
the external source, the expectation value evolves in time as a solution of the effective equations of motion in the 
medium. For a detailed description of this method in non-equilibrium quantum field theory see refs. |26ll45ll46l| . The 
external sources rja in the Lagrangian density (|II.ip induce an expectation value of the neutrino field 

■ipa = {va) =T'rpiya (11.13) 

where p is the full density matrix of the medium. In linear response this expectation value is linear in the external 
source and obeys the effective Dirac equation of motion in the medium ;4(l]. It is most conveniently written in terms 
of the spatial Fourier transforms of the fields, sources and self-energies ipa{k, t),r]a{k, t), I](fc, t — t') respectively^. The 
one loop self-energies with neutral and charged current contributions had been obtained in refs.[22ll25ll26j|. and the 
effective Dirac equation in the medium up to one loop has been obtained in the real time formulation in ref.|2g. It 
is given by 



*7°^ - 7 • fc ) -^afc - Mab + Kl^L 



-t 
Mk,t)+ I dt'j:ab{k,t-t')LMk,t')^-f]a{k,t), (11.14) 



where L is the projector on left handed states, S^'^'^L is the (local) tadpole contribution from the neutral currents and 
Y,ab{k,t — t') is the spatial Fourier transform of the (retarded) self energy which includes both neutral and charged 
current interactions, and whose spectral representation is given by 

/OO TT 
— ^ImS(fc.fco)e-*'''«(*-*') ; I](fc, fco) = SH'(fc, fco) + Sz(fc. fco) (11-15) 



^ We have kept the same functions to avoid cluttering the notation, but the label k makes it clear that these are the spatial Fourier 
transforms. 



where we have separated the charged and neutral current contributions respectively. 

The external source term 77 allows to "prepare" a determined initial state, leading to the time evolution of ip a-s an 
initial value problem. This approach is implemented as follows. Consider switching on the source adiabatically from 
t = — cxD up to i = and switching it off at t = 0, 



Vaik,t)^r]a{k,0)e''e{-t) ; e ^ 0^ 



(11.16) 



ft is straightforward to confirm that for < < the solution of the Dirac eqn. pi.l4|l is given by 



Mk,t<0)=Mk,0)e' 



(11.17) 



Inserting this ansatz into the equation (|II.14|I yields a linear relation which determines the initial value 4'a{k,0) 
from r]a{k, 0), or equivalently, for a given initial value ip{k, 0) allows to determine the adiabatic source that yields this 
initial value problem. The evolution for i > is determined by the following effective (retarded) Dirac equation, 



i7"^-7-fc) Sab-Mab + Kl^L 



Mk,t)+ / dt'Y.ab{k,t~t')Li}jb[k,t') = -Va(fc,0) 



dko Iml](fc, ko) 

TT fco 



'ikot 



This equation can be solved by Laplace transform. Introduce the Laplace transforms 



(11.18) 



tpb{k,s) 



'Mk,t) ; I](fc,s) 



^E(fc,t)= ^ 



dkci lmS{k, fco) 



TT ko ~ is 



(11.19) 



where we have used eqn. (|II.15ll to obtain the Laplace transform of the self-energy. At this stage it is convenient to 
establish contact with the more familiar description of the Dirac equation in the frequency representation (Fourier 
transform in time) by introducing the time Fourier transform of the retarded self energy, 



S(fc,w) 



dko ImE(fc, fco) 



TT kf) — uj — ie 
related via analyticity to the Laplace transform, namely 

E(fc, s) = E(fc, u) — is — ie) 
In terms of Laplace transforms the equation of motion becomes the following algebraic equation 



(11.20) 



(11.21) 



Dab{k,s)ipb{k,s) =i [j^ 6ab + — ^ab{k,s) -Y.ab{k,0) L]ipb{k,0) 

IS 



(11.22) 



where D{k,s) = D{k,uj — is — ie) is the analytic continuation of the Dirac operator in frequency and momentum 
space 



Dab{k, UJ) = [(^fu ^l-k^Sab- Mab + Kt L + S,b(fc, io) 

The full space-time evolution of an initial state is determined by 



L 



(11.23) 



ds 



M^,t)^ /d^fce*^^ / —Djik,s)itj')Mk,0)e'\ 



2m 



(11.24) 



where F is the Bromwich contour in the complex s plane running parallel to the imaginary axis to the right of all the 
singularities of the function '0(fc, s) and closing on a large semicircle to the left. We have simplified the expression for 
the eqn. (|II.24|) by discarding a perturbatively small correction to the amplitude of 0{Gf), given by the self-energy 
corrections on the right hand side of eqn. (|II.22|I . Therefore the space-time evolution is completely determined by 
the singularities of the function '(/'(fc, s) in the complex s-plane. Up to one loop order and for temperatures well below 



the mass of the vector bosons, the only singularities are simple poles along the imaginary axis, corresponding to the 
dispersion relations of the propagating modes in the medium. In this temperature range absorptive processes emerge 
at the two loop level, consequently these are of 0{Gp) and are neglected in the one loop analysis presented here. The 
integral along the Bromwich contour in the complex s-plane can now be written 



7 /"OO 7 



2ni 



27r 



(11.25) 



where the frequency integral is performed along a line parallel to but slightly below the real lu axis closing counter- 
clockwise in the upper half plane. 

The one loop contributions to the self-energy for w, fc, T <C Mw were obtained in reference |23. l25l l2a| and found to 
be of the form,26j 



Elf + Eafc(fc, c^) = 7°A(co) - 7 • kB(fc) 



(11.26) 



where A(w) and B(fc) are 2x2 diagonal matrices in the neutrino flavor basis. ^. To lowest order in k/Mw ;uj/Mw 
these matrices are found to be|ilO, 22, 25 , 26] 



A(w) 



Ae{uj) 

Af^iiu) 



B(fc) 



B^{k) 
B^{k) 



(11.27) 



Imposing charge neutrality, the results of reference |2g (see also[T3 12^ ^^ Ull) are summarized in the following 
regimes: high temperature and density, relevant during the early Universe or low temperature and high density, 
relevant for cold dense matter in supernovae, neutron stars or the sun 

• me < T < m^ 

In this regime we consider the following degrees of freedom v, e and p, n in nuclear statistical equilibrium. The 
matrix elements are given by 



Aeiuj) 

Be{k) 
B,(k) 



Gp n^ 

Gp n-y 



c. 



Ttt* ljT 



C 



60C(3) M^ 



^60 C(3) M, 



(2 -I- cos^ 



cos 9,, 



Gp n^ 7 n^ kT , ^ ^ 

Gp rij 1 -K^ kT 2 a 



V2 180 C(3) M, 



w 



where in terms of the asymmetries Lf with / = VjejP^n 



(11.28) 

(11.29) 

(11.30) 
(11.31) 



i-'e — ^ -t^i/e ^^ -^1^^ ~ 3 -Le ~ J^n ', ~ 1-^ — T^^i^ii. '^ ^t^c ~ ^" 



(11.32) 



• me, m^ < T < Mw 

In this regime the relevant degrees of freedom are the lightest deconfined quarks u, d and leptons 



The equivalence with the notation of ref. '22!] is (see eqn. (2) in ref . [2^ ) : ajvH = B(fe)/fc ; fejvij = A(aj) — lo B(A;)/fc. 



Ae{0j) = 


Gpn^ 
V2 






M^) = 


Gf n^ 

V2 


~ Ttt-* wT , 2. N 




Be{k) = 


Gfu^ 77r4 fcT / 2/, 60 /me\2 
V2 180C(3) Ml y ' ^°" "^ I^^Kt) 


B^.{k) = 


Gfu 

V2 


. 77r4 fcT / 2., 60 





(11.33) 
(11.34) 

(11.35) 
(11.36) 



where 



L^^ + U^ - 3Le + (1 - 4sin2 B^){2L^ - L^,) - (1 - Ssin^ 9^)Lu - 2i, 



M - --L^„ +L^^ -3i^ + (l-4sin"0^)(2i, -L^)-(l-8sin^6'„)L„-2Ld 



(11.37) 

(11.38) 



• Cold dense matter with e, v, ^ItPtU : 






GpUe 

V2 



A„ 



GfM^ 
V2 



ACe+AA,^-3AAe-AA„ 
K^ + K^ - ■N'n 



Be,M = 



(11.39) 
(11.40) 
(11.41) 



• Cold dense matter with quarks and leptons : 



Ae 
-Me 






V2 
-^AC, + AA,^ - 3JVe + (1 - 4sin2 0^)(2AAe - A^^) - (1 - 

-^M,^ + AC. - 3M^ + (1 - 4sin2 e^){2Me - M^) - (1 



8sm^0^)Mu-2Afd 
-8sm^e^)Afu-2Afd 



(11.42) 

(11.43) 
(11.44) 



where n^ = 2 C(3) T^/ir'^ is the photon density, and 



Lf = 



Mf 



' ^f = ^I^ANfip)-NM 



Aff are the particle-antiparticle asymmetry densities. 



(11.45) 



B. Dispersion Relations and Mixing Angles in the Medium 



The simple poles of the integrand in ljll.25|l are the solutions of the homogeneous Dirac equation 

7° wl - 7 • k fcl - M + (i"A{lj) - 7 • kB{k)] L ^/j{cu, fc) = , 



(11.46) 



where 1 is the 2x2 identity matrix in the flavor basis in which the fleld ip{uj,k) is given by 



ip{uj,k) 



_ ( Vei'^-.k) 

i'f,{uj,k) 



(11.47) 



with Ve{i-^, k) and Vfj,{u), k) each being a 4-component Dirac spinor. 

It turns out to be most convenient to work in the chiral basis in which the left-handed and right-handed components 
of the Dirac doublets are written as M 



*'=!:„ 



h=±l 







; ^R=Y. 



h=±l 







(II.. 



where the two component Weyl spinors v^^'^ are the eigenstates of the helicity operator a ■ k with eigenvalues h = ±1. 
To the leading order in Gf, the left-handed flavor doublet 



ip^^Kuj,k) 



vi''\LO,k) 
4''\LU,k) 



(11.49) 



obeys the following effective Dirac equation |26 



[{uj^ - fc2)l + {uj- hk){A + hB) - M^] ip^''\uj, k) = 
while the right-handed doublet is determined by the relation |26l | 



(11.50) 



? 



CO 



{^,k) = -^^^^±^M^('^\u,k) 



fc2 



(11.51) 



The propagating modes in the medium are found by diagonalization of the above effective Dirac equation. This 
can be done by performing a unitary transformation ip^^'{uj, k) — Um X {'^^ ^) where 



Uih) 



cosBt^ sin^^'*) 



W ...flW ) ' X^'n^,fc)= 1 ,f/.) 



sin Qm cos 



_ { vf\uj,k) 



(11.52) 



and a similar transformation for the right handed doublet ^'''•'(cij, /c), with the medium mixing angle Om depending 
on h, k and to. Upon diagonalization, the eigenvalue equation is given by |2g| 



L2 - fc2 + i Suico, fc) - i {Ml + Ml) -ISM^ n^ico, k) [l \^ W"'(^, k) - 
where Shito, k), /^.^{lo, k) and ^^{lo, k) are respectively given by 



(11.53) 



Sh{uj,k) = {u-hk)[A,{uj)+A^{uj) + hB,{k) + hB^{k)] , 

A;,.(cu,fc) = {uj-hk)[A,{uj)-A^{oj) + hB,{k)-hB^,{k)] , 

1 
2 



nh{uj,k) = 



cos 29 



Ahiuj,k) 
SM^ 



sin^ 26 



This requires the matrix elements in Um to be of the following form 



(h) 



(11.54) 
(11.55) 

(11.56) 



. „^n,w ,N sin 20 

sm2e!^>{uj,k) = TT-. — 7T ; cos 



^h(w,fc) 



„„„9/l Afc(aj fc) 

\lh(uj,k) 



(11.57) 
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A resonance is available whenever 



»..-M^,.o 



(11.58) 



in which case sin 26'm(w,fc) ~ 1 ; cos29m {uj,k) ^ 0. 

The solutions of eqn. (|11.53|l yield the dispersion relations uj{k) of the "exact" quasiparticle states in the medium 

and correspond to the "exact poles" of the Dirac propagator. The dispersion relations wi (fc^X) for the propagating 
modes in the medium are found in perturbation theory consistently up to 0{Gf) by writing |26| 



where Ea{k) = ^k? + M^ , and SuJa {k, A) are found to be 



a= 1,2 ; A = ± 



(11.59) 



5c.i'')(fc,A) = --^lsh{XEa{k),k) + {-lySM^ {nh{\Ea{k),k) - 1) 



(11.60) 



For relativistic neutrinos with k ^ Ma the dispersion relations uJa{k) ; a = 1,2 for the different cases are given 
to leading order in Gp by 

• Positive energy, negative helicity neutrinos, A = +l,h — —1: 



-(^)-^+f -i 



S^ik, k) + {-lY5M'^Ul_{k, k) 



• Positive energy, positive helicity neutrinos, A = +1, h = +1: 



UJa{k) = k + 



Ml_l_ 
2k Ak 



S- 



-{k, k) + {-IfSM"^ (p.+ {k, k) - l) 



Lo — hk 



M 
'2k 



where we have neglected corrections of order SM"^ / M . 

Negative energy, negative helicity neutrinos, A — — l,/i = 1: 



Af^ 1 



Oa{k)^-k--^ + 



2k 4fc 



S- 



{~k,k) + {^iydM^(n^{-k,k) - 1) 



hk: 



M 
'2k 



where we have again neglected corrections of order SM'^ / M . 
• Negative energy, positive helicity neutrinos, A = —l,h — +1: 



(11.61) 



(11.62) 



(11.63) 



Mi . 1 



UJa{k)^-k~-^^ 



2k 4fc 



5+(-fc, k) + {-lySM^ h+i-k, k) - l\ 



In the above expressions the Vt± are given by eqn. (|II.56p . 

The vacuum and medium oscillation time scales are respectively defined as 



(11.64) 



T — 



27r 



E^ — Eo 



-^ Tiled — 



2tt 



Xh) 



c^r(fc,A)-4"^(fc,A) 



(h). 



(11.65) 



In the relativistic case when k ;:?> Ma, we find 
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Ank 



A-Kk 



5M'^Vth{\k,k) 



(11.66) 



leading to the relation 



Tmed _ sin26'm (c^,fc) 
T^ac sin 20 



(11.67) 



III. SPACE-TIME PROPAGATION OF A NEUTRINO WAVE-PACKET. 

We now have all the ingredients necessary to study the space-time evolution of a initial wave packet by carrying 
out the integrals in eqn. l|II.24(l . For this purpose it is convenient to write ^{k, 0) = i^R{k, 0) + i^L{k, 0) and expand 
the right and left-handed components in the hclicity basis as in eqn. (|II.48p . namely 



Mk,0)= J2 



h=±l 



u('')®^('')(fc,0) 



; ^fl(fc,0)= Y. 



h=±l 







(IILl) 



where 



^('')(fc,0) 






(fc,0) 
(fc,0) 



C''^(fc,o) = 






ih) 



(fc,0) 

(fc,0) 



(III.2) 



The general initial value problem requires to furnish initial conditions for the four components above. However, 
an inhomogeneous neutrino state is produced by a weak interaction vertex, which produces left handed neutrinos, 
suggesting to set v^^^ (k, 0) = 0; v ^ (k, 0) = 0. Without loss of generality let us consider an initial state describing an 

inhomogeneous wave-packet of electron neutrinos of arbitrary helicity, thus i^g^ (fc, 0) 7^ 0; vl^lik, 0) = 0. 

In the cases of interest neutrinos are relativistic with typical momenta k 3> M . Following the real time analysis 
described in detail in ref.l26l in the relativistic case we find 



(h), 



^'''>ik,t) = -:.ll'ik,0) 



l+C 



(h) 



'S_ 



(h) 

-h 



,-iuj'-l'\k-h)t 



+ 



l-C 



ih) 



s 



ih) 

-h 



^-^J^\k,-h)t^Q(M_\ 



(IIL3) 



^('')(fc,t) 



,{h) 



'eh 



{k,0) 




1 + cos 20 
- sin 29 



^-luj["\k,~h)t ^ 



'iJ^\k,h)t 



1- cos 26* 
sin 26* 






V k 



(IIL4) 



where (p^^^'lkji) and £,^^^{k,t) are the flavor doublets corresponding to the left-handed and right-handed neutrinos 
with helicity h respectively. The upper component corresponds to the electron neutrino vi {k, t) while the lower 



component corresponds to the muon neutrinos i^jj, \k,t). The factors Cj^ (fc) and 5j^ (fc) are defined as 



{h). 



:('')/ 



cf)(fc)=cos[20(^)(Afc)l ; 5f'(fc)=sinte)(Afc)' 



(IIL5) 



The suppression factor M/k in the right handed component ljIII.4|) is of course a consequence of the chirality fiip 
transition from a mass term in the relativistic limit. For relativistic neutrinos and more specifically for neutrinos in 
the medium prior to BBN with fc ~ T ~ few MeV the right handed component is negligible as expected. 
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The one-loop computation of the self-energy performed above does not include absorptive processes such as collisions 
of neutrinos with Icptons (or hadrons) in the medium. Such absorptive part will emerge in a two loops calculation and 
is of 0{G%). While we have not calculated these contributions it is clear from the analysis what it should be expected: 
the frequencies ciJi_2(fc) are the "exact" dispersion relations of the single particle poles of the Dirac propagator in the 
medium. At two loops the self energy will feature an imaginary part with support on the mass shell of these single 
particle states. The imaginary part of the self-energy evaluated at these single particle energies yield the width of 
the single quasi-particle states ri(A:);r2(fc) and the oscillatory exponentials in the expressions above are replaced as 
follows 



-iua (,k) t 



-ra{k)t -iUa{k)t 



a= 1,2 



(III.6) 



While our one loop calculation does not include the damping rates Ta we will invoke results available in the 
literature[i3|23i|2J,|23 to estimate the coUisional relaxation time scales (see section Hvjl . 

The corresponding fields for the left-handed and right-handed component neutrinos in configuration space are 
obtained by performing the spatial Fourier transform 



V>"'\r,t) 



(2^ 



(^('')(A?,i) e 



ik-r 



(III.7) 



&\r,t) 



(fik 



eW(fc,t)e' 



ik-r 



(III.8 



For an arbitrary initial configuration these integrals must be done numerically, but analytic progress can be made 
by assuming an initial Gaussian profile, describing a wave-packet in momentum space centered at a given momentum, 
fco with a width a. While the width could generally depend on helicity we will consider the simpler case in which it 
is the same for both helicities. Namely, we consider 



.ih)i 



.(h), 



^eLik,0)-'^'el'iO)[^ 



exp 



(fc - fco)^ 

4cr2 



(III.9) 



where z/^^ (0) is an arbitrary amplitude and assume that wave packet is narrow in the sense that a <C fco . In the limit 

CT ^ the above wave-packet becomes I'g^ (0)(5^(fc — fco). In the opposite limit of large a the wave packet describes 
an inhomogeneous distribution spatially localized within a distance ~ l/a. For a narrow wave packet the momentum 
integral can be carried out by expanding the integrand in a series expansion around fco keeping up to quadratic terms. 



The typical integrals are of the form 



A. Integrals 



/(f,i) 



TT \ 2 



_fk_ 

(27r)3 



^(fc) exp 



(fc~fco)^ 

4cr2 



-|- ik ■ r — iLu{k) t 



(III.IO) 



where A stands for the factors (1 ± C) ; 5 in eqns. (jIII.3IIII.4|) . and uj{k) are the general dispersion relations obtained 
above. The computation of these integrals is simplified by noticing that for any function F{k) the expansion around 

fco up to quadratic order is given by 



F{k) = F(fco) + i^'(fco) ko • (k - ko) + i (F"{ko) Pll (ko) + ^^ P.^(ko) ] (k - ko).(k - ko), 



(III.ll) 



where 



Fl(k)-k,k, ; i^j(k)-<5,,-k,k, 



(III.12) 
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The result of the integration can be written more compactly by introducing the following quantities 



aht) = a^\ ^^ci,||(t) 1-4 



(2 



(III.13) 



al{t) = a^ \ '^^} ^<^^[t)(l-^- 



t2 



T± 



where we have introduced the perpendicular and parallel dispersion time scales given respectively by 



T± 



fco 



1 



2a^Vo 



ni 



2cr2tj"(fco) 



7^T_l 



(III.14) 



It will be seen in detail below that these two time scales are indeed associated with the spreading of the wave packet 
in the transverse and longitudinal directions. 

The group velocity Vg and effective Lorentz factor^ 7 are given by 



Ug=w'(/co)ko ; 7^ 



/cow"(fco) 
The transverse and longitudinal coordinates are 

-^11 (0 = ko [r- ko - Vgtj ; X^^r-'ko^r- koj 

and in terms of these variables we finally find 



(III.15) 



(III. 16) 



I{f,t) = 



a\\{t)al{t) 



A{ko;r,t) e<feo-r--*(r,t)t) g-(<i>_L(t)xi+<E>||(t)x^(t)) 



where the phase 



(III.17) 



^{f,t)^iu{kQ) + 



<&i(i) ^2 , *llW 



X 



xfAt) 



(III. 18) 



and 



A'iko) 



A{ko;r,t) = A{ko) + 2tA'{ko)aj{t)ko ■ X\\{t) + A"{ko)al{t) (l - 2alit)Xl) + ^^^^<Ji\it) (l - 2a|(i)X|f (t)) 

(III.19) 
Neglecting the prefactor y^(fco; J^, i) we see that 



|/(r,t)p« 



('4)('4r 



~2(*i(t)Xi+*l|(t)Jf2(t)) 



(III.20) 



describes a wave-packet moving in the direction parallel to the momentum feg with the group velocity Vg and dispersing 
both along the perpendicular and parallel directions. The expressions for ^±{t) and ^\\{t) given by eqn. IJIII.13|) 
clearly show that the dispersion time scales along the parallel direction and transverse directions are given by T|| , t± 
respectively and T|| displays the time dilation factor 7. The wave packet is localized in space within a distance of 



For the usual dispersion relation a;(A;) = \/k^ + M^ it is straightforward to confirm that 7^ = (1 — v'^) ^ 
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order l/-\/$(i) oc 1/ \/a in either direction. Small a localizes the wave packet in momentum space while large a the 
wave packet is spatially localized. For large a the integrals must necessarily be performed numerically. 

This discussion highlights that the derivative terms in the prefactor ^(fco;r, i), which are a consequence of the 
momentum dependence of the mixing angles correspond to an expansion in the ratio cr/fco. This can be understood 
from the following argument: A ^ (1 ± C), 5, hence its derivatives with respect to momentum are of the form f{k)A' 
with f{k) being smooth and bounded functions of 0(1), while A is at most of the form Aofc + Aifc^ in the relativistic 
limit, (see eqn. (|II.55(I ^ therefore A' m A/k. These derivatives multiply powers of cr_L,||^±.||, and the exponential 
damping in / restricts these contributions to the range |crj_ ||Xj_ ii | « 1. Therefore in the narrow packet approximation 
cr <C fco the higher order derivative terms are suppressed by powers of cr/fco ^ 1- We have invoked this narrow packet 
approximation to carry out the momentum integral, therefore consistently with this approximation we will only keep 
the first derivative term, which is of ©(a/fco) and neglect the higher order derivatives, which are of higher order in 
this ratio. Namely in the analysis that follows we approximate 



\Aiko;r,t)\^^\A{ko)\' 



1 + 4 



^^(fco) 
A{k„) 



<D||(t)ko-X||(i) 



(III.21) 



In this manner we consistently keep the lowest order corrections arising from the momentum dependence of the 
m,ixing angles in the medium. 

We now have all the ingredients for our analysis of the space time evolution. The above general expressions for the 
time evolution of initially prepared wave-packets, eqns. piI.3IIII.4ll combined with the dispersion relations obtained 
in section l|IIB|) provide a solution to the most general case. We focus our discussion on the case of the early 
Universe, in which the typical neutrino energies are ~ MeV. With (active) neutrino masses in the range Ma ~ eV and 
SM"^ ^ 10"'°^ — 10^'^ it is clear from the results above that the amplitude of the right handed component is suppressed 
by a factor M/fc r^ 10~^ and the medium corrections to the dispersion relations for positive energy neutrinos with 

positive helicity and negative energy neutrinos with negative helicity are suppressed by a factor M /fc^ with respect 
to the opposite helicity assignement. Therefore in what follows we focus our discussion to the case of left handed 
negative helicity neutrinos (and positive helicity antineutrinos). 



B. Space-time evolution and oscillations. 

We now focus on describing the evolution of negative helicity neutrinos or positive helicity antineutrinos. 

The initial state considered above corresponds to a wave-packet of electron neutrinos at t = but no muon neutrinos. 
The lower component of the flavor spinor in eqn. IJIII.3IIII.7|I describes the wave-packet of the muon neutrino at any 
arbitrary time. We begin by studying the transition probability from an initial electron neutrino wave packet of 
negative helicity to a muon neutrino wave packet. 

Using the results obtained in the previous section for the integrals in the narrow packet approximation we find the 
transition probability 



.W, 



5(fco) rii 



Pe^,(f,i) = i<i:(f,i)r = jW':i>{k,or\siko)[ 

|/i(r,t)|2 + |/2(f,t)|2 - 2 |/i(f,t)| |/2(f,t)| cos[(vI/i(r,t) - *2(f,t)) t] 



(III.22) 



where S — sin 



26'm'(±fc) and /i,2(^, ; ^1,2 correspond to the integrals and phases given by eqn. Ijlll. 171111. 18() 

with the frequencies 0Ji,2{k) for negative helicity given by eqns. (|II.61|I . In the expression above we have taken a 
common prefactor by neglecting the differences between the group velocities and the masses, taking Vg — I, and 
^jl,?^!! correspond to *I'||,T|| with a mass M. We focus our attention on the interference term which is the space- 
time manifestation of the oscillation phenomenon and features the oscillatory cosine function. The amplitude of the 
oscillation, |/i(r, t)l2{r, t)\ describes the product of two wavepackets of the form given by eqn. (|111.20p . 
It is convenient to write the product |/i /2I in the following form 



\h{r,t) h{r,t)\ 






-(*i,l(t) + *i,2(t))^i p-'i'CM(t)Xl_.^jit) „-*H(t)Xj(t) 



(III.23) 
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where we have introduced the center of mass (CM) and relative (R) variables 

XcM = kQ[r-ko-vcM[t)t] ; vcAiit) = — — -r (III. 24) 



$11 1(0 + ^112(0 

Xr = X|| 1 - i|| 2 = - (Vg 1 - Wg 2) i (in.25) 

$CM = $11 1 + $11 2 ; $i? = ^ " , r (III.26) 

*P|| 1 + *P|| 2 

The integral ljIII.23p describes the product of two gaussian wave packets spreading in the transverse and longitudinal 
directions and separating in the longitudinal direction because of the difference in group velocities, made explicit by 
the term $ij(i)X2 (t). 

The first two terms in eqn. (|III.22|I describe the incoherent sum of the probabilities associated with separated wave 
packets of propagating mode eigenstates , in the third, interference term, the product I/1II/2I is the overlap between 
these two wave-packets that are slowly separating because of different group velocities. As discussed above a two loop 
calculation of the self-energies will lead to a quasiparticle width and a damping rate Ta for the individual quasiparticle 
modes of frequency 0Ja{k), the discussion leading up to eqn. (|III.6|I suggests that the integrals 



\Iair,t)\^e-^-'^'^^'\Ia{f,t)\. (III.27) 

1. Coherence and "freeze-out" 

Since Xr — {vg2 ~ Vgijt does not depend on position, the overlap between the separating wave packets becomes 
vanishingly small for t » tcoh where the coherence time scale tcoh is defined by 

'^B{tcoh){Vg2~Vg^ftl^^ - 1 (III.28) 

Before we engage in an analysis of the different cases, it is important to recognize that there are several dimensionless 
small ratios: i) cr/fco ^ 1 describes narrow wave-packets, this approximation was implemented in the calculation of 
the integrals, ii) M /k <^ 1 in the relativistic limit with k ^ MeV for example in the early Universe near the epoch of 

BBN or for supernovae, iii) 5M'^ / M <^ 1 describes a nearly degenerate hierarchy of neutrino masses. Since in the 
relativistic limit fig — f2g ^ SIvP /k"^ we can neglect the difference in the masses in ^y and write $|| 1 ~ $|| 2 ~ *&|| 
where the masses are replaced by the mean mass M given by eqn. (|II.9|I . and similarly for <i>^. Therefore to leading 
order in small quantities we can replace $jj above by $11/2 leading to 

— 2cr^ 1 

$CAf(i) = 2$||(t) = - jT ; vcM^7^{vgi+Vg2) (111.29) 



where t\\ is given by eqn. (|III.14|I for M, and 



1- 



(i) 



-$ll(t)X^(0= ^7,2 (III-30) 



2 



1 



where we have introduced the time scale 



(T\Vg2~Vgl\ 

The coherence time scale is the solution of the equation 



\/2 
tc = -, r (111.31) 



Ml-. 



_ 2 



(III.32) 
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The expression (jIII.3011 reveals a remarkable feature: for t ;:!> r|| the overlap between the separating wave-packets 
saturates to a time independent value 



-$||(i)X^(t)^(^^j . (III.33) 

This effect has been recognized in ref. '37!| and results from the longitudinal dispersion catching up with the separation 
of the wave packets. This phenomenon is relevant only in the case when tc > T|| in which case the overlap of the 
separating wave packets "freezes" and the packets maintain coherence for the remainder of their evolution. There are 
two distinct possibilities: 



In case (a) we can approximate 



tc « r|| : (a) (III.34) 

tc » r|| : (b) (III.35) 



i<i>ll(t)X]-.(t)«('A) (III.36) 



since during the time interval in which the separating packets maintain coherence i <C ic ^ t\i and in this case the 
relevant coherence time scale is tc- 

In case (b) the "freeze out" of coherence results and the long time limit of the overlap between the wave packets 
in the longitudinal direction remains large and determined by eqn. Ijlll.33|) . 

However, while this "freezing of coherence" phenomenon in the longitudinal direction ensues in this regime, by the 
time when the coherence freezes t r^ t|| the wave packet has spread dramatically in the transverse direction. This is 
because of the enormous Lorentz time dilation factor in the longitudinal direction which ensures that t ^ tu 3> tj_ (see 
eqn. (|111.14p 'l. The large spreading in the transverse direction entails a large suppression of the transition probability 



r\\ J 7^ V fco 



.(r,i^r||)cx^ --4-(-) ■ (111-37) 



For M ~ eV and fco ~ MeV the above ratio is negligible. Therefore while the phenomenon of freezing of coherence 
is remarkable and fundamentally interesting, it may not lead to important consequences because the transition prob- 
ability is strongly suppressed in this regime. Therefore in the time scale during which the transition probability is 
non-negligible, namely t <C tii the overlap integral can be simplified to 



4*li(t)xa(t) ^p-ii:)" 



2. Effective oscillation frequency 



(III.38) 



Another aspect of the interference term is the effective time dependent oscillation frequency ^i(r, t) — ^2(^1^) 
where the ^a are given by eqn. (|111.18|l for the frequencies uja{k) of the propagating modes. The spatio-temporal 
dependence of this effective phase is a consequence of the dispersion of the inhomogeneous configurations, encoded in 
the functions $ and results in a drift of the oscillation frequency, a result that confirms a similar finding in the vacuum 
case in ref. j49j . Because of the exponential fall off of the amplitudes the maximum value of the drift contribution 
is achieved for $_L,a^i ~ 1 ; ^||,o-'^|?a(^) ~ ^' namely in front and back of the center of the wave-packets, both 
in the transverse and the longitudinal directions. Furthermore, because of the Lorentz dilation factor, T|| 3> t± for 
relativistic neutrinos. Therefore we can approximate the effective oscillation frequencies as 

2cr2 
*i - *2 - LOlik) - U2{k) -f ^— (%i - Vga) (III.39) 

Kg 

The dispersion relations and mixing angles obtained above along with the the results (|111. 31111. 4ll . yield the complete 
space-time evolution for wave-packets with initial conditions corresponding to an electron neutrino. Rather than 
studying the general case, we focus on three different situations which summarize the most general cases, i) A;j/i5M^ <C 
1 corresponding to the case of vacuum oscillations, ii) A^/SM'^ ~ cos 20 corresponding to a resonance in the medium, 
and iii) Ah/SM'^ 3> 1 corresponding to the case a hot and or dense medium in which oscillations are suppressed. 
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C. ^h/5M'^ <^ 1: vacuum oscillations. 



We study this case not only to compare to results available in the literature, but also to establish a "benchmark" 
to compare with the results with medium modifications. Beuthe pT\ has studied the propagation of neutrino wave- 
packets in the vacuum case including dispersion and in ref. [43 an effective frequency similar to (jfff.l8|) has been found 
for wave-packets propagating in the vacuum. In this case for positive energy, negative helicity neutrinos with a ~ 1,2 



Wa(fco) ~ fco + 1^ ; Vg,a 



2ko 



Ml 
'4 



1 a "/ 7 \ a 



ap 

.3 




leading to the vacuum time scales 



^c,v 



2V2k^ 
a\SAP\ 

jL3 

2am^ 



In the case when 



the vacuum coherence time is given by 



ko SJVP 



AM 



>1 



^c,v — "^11 



a AM 



fco (5M2 



2Kq 



a\5M'^ 



<T,, 



(III.40) 

(III.41) 
(III.42) 

(III.43) 
(III.44) 



and the overlap between the separating wave-packets vanishes well before the packets disperse appreciably along the 
longitudinal direction. On the other hand, in the case when 



kQ SM^ 



AM 



<1 



(III.45) 



the spreading of the wave packets catches up with the separation and the overlap between them freezes when t — 

2 2 

tf^y — T|| ~ k'^/2a^M . With \5M'^\/AM ^ lO"'* for solar or ~ lO^'^ for atmospheric neutrinos the phenomenon of 

"freezing" of the overlap and the survival of coherence is available for 



fco 4 M 



(III.46) 



which is well within the "narrow wave packet" regime. However, as discussed above, when the coherence freezes the 
transition probability has been strongly suppressed by transverse dispersion. Therefore during the time scale during 
which the transition probability is non-negligible we can approximate the exponent in ()III.23|) 



The effective oscillation frequency is given by eqn. Ijlll.39|) which becomes 



(III.47) 



5M^ ( 2^2 
*l-*2 r^ I 1- 



(III.48) 



while the corrections tend to diminish the oscillation frequency, these are rather small in the narrow packet approxi- 
mation. 
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D. Medium effects: near resonance 



In refs . [10|, l22l |24L l25l l26l | it was established that if the Icpton asymmetries are of the order of the baryon asymmetry 
r] ^ 10~^ there is the possibiUty of a resonance for the temperature range rrig <C T ^ m^ for positive energy 

2 

negative hehcity neutrinos with ui{k) ^ k + AI /2k ; ft, = —1 or positive energy positive hehcity antineutrinos with 
Lu{k) ^ —k — M /2k ; h = \ respectively. It is convenient to introduce the following notation 



^5 = 






eV 



(III.49) 
(III.50) 



// the lepton and neutrino asymmetries are of the same order of the baryon asymmetry, then 0.2 < |£g| < 0.7 and the 
fitting from solar and KamLAND data suggests jf^sj ~ 8. In this temperature regime we find|2g for positive energy, 
negative helicity neutrinos 



A_(fe,fc) 



4 



O.IT 

MeV 



k 
T 



-Co 



( 2T 



VMeV 



(III.51) 



and for positive energy positive helicity antineutrinos 



A+{-~k,k) _ 4 /O.lTy k 
SM^ ~ ^ \Mev) T 



Co 



2T \ k 
MeV J T 



(III.52) 



In the above expressions we have neglected terms of order -^. With k ^ T and in the temperature regime just 
prior to BBN with T ~ few MeV the lepton asymmetry contribution C is much smaller than the momentum dependent 
contribution and will be neglected in the analysis that follows, therefore we refer to A^(Afc, k) and Shi^k, k) as A(fc) 
and S{k) respectively since these are independent of h, A in this regime. In this temperature regime we find for both 
cases (negative helicity neutrinos and positive helicity antineutrinos) the following simple expressions 



A{k) 



567r^ Gpk^T'^ 

45n/2 



M^ 



S'(fc) w A(fc)(l + cos2 6i^). 



(III.53) 



A resonance is available when A(fco) ~ (5M^cos20, which may occur in this temperature regime for fcg '^ ^^ ~ 
3.6 MeV 10, 22, 24, 25, 26] for large mixing angle (^12) or A; ~ T - 7 MeV for small mixing angle (6*13). Near the 
resonance the in-medium dispersion relations and group velocities are given by 



Wa(fc) 



"g,a 



fc + — 2- — <^ (1 + cos2 e^) COS 29 + (-l)'^-i (1 - sin 261) 

2k 4fc I 

/V/f2 r»,r2 ( 

1- ^ + ^^ ai +008^0^) COS 20+ (-1)^-1 (1-sin 20) 



(III.54) 
(III.55) 



Again we focus our discussion on the interference terms in the transition probability (|III.22|I . in particular the 
medium modifications to the oscillation frequencies and coherence time scales. To assess these we note the following 
(primes stand for derivatives with respect to k): 



n{k)\res ^ Sm2e ■,n{k)\res^O ; ^ (fc)|r 



4 cos^^ 
fc2 sin 29 



(III.56) 



which when combined with equation (|II.61|I yield 



Wi(fc) — UJ2{k) 



SAP 
~2k 



■ sin 29 ; v„ i — v, 



9,2 



SAP 



■sin 20 



(III.57) 



We also note that near the resonance 



sin' 29„i{k) ex cos20„ 



0, 
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(III.58) 



therefore the corrections arising from the energy dependence of the mixing angle in the transition probability pil.22|l 
become vanishingly small. The transverse and longitudinal dispersion time scales are given by 



fcn 



2cr^ V„n 



Ta\\ 



l-(-l) 



a-iSAP l+COs2 2e' 



2M 



sin 26* 



II 



Kg 



2<t2m 



(III.59) 



Therefore in the medium near the resonance, the argument of the exponential that measures the overlap between 
the separating wave packets is given by 



-$11 {t)Xl{t) ^ -^ ^ sm 20 =j-- 



where 






T^ctn 



2Kq 



(T|(5M2|sin26l sin 26* 
The effective oscillation frequency (jIII.39|) is given by 



(III.60) 



(III.61) 



*1 -*2 



5M^ . „ / 2(72 
— — sin 26* 1 - —2- 



(III.62) 



which when compared to the vacuum result ljIII.48|) confirms the relation between the vacuum and in-medium oscil- 
lation time scales (lii.67|l since near the resonance sin26'„ ~ 1. 

We conclude that the main effects from the medium near the resonance are an increase in the coherence and in 
the oscillation time scale T = 2tt/\'^i — ^2! by a factor l/sin26'. For solar neutrino mixing with sin26'i2 ~ 0.9 the 
increase in these time scales is at best a 10% effect, but it becomes much more pronounced in the case of atmospheric 
neutrino mixing since sin 2^13 <C 1 . 



1. Ah/5M^ ^ 1 oscillation suppression by the medium. 



In the temperature or momentum regime for which Ah/SAI'^ ^ 1 the expression for the in-medium mixing angles 
(|II.57|I reveals that cos26'm ^ — 1. In this case the in- medium mixing angle reaches 9m — > 7r/2 and the transition 
probability Ve^fj, vanishes. Eqn. IjIII.Sp shows that in this case an electron neutrino wavepacket of negative helicity 
propagates as an eigenstate of the effective Dirac Hamiltonian in the medium with a dispersion relation 



uj2{k) ^ vk+ —^ 



1- 



45 a/2 M^ 



(1 -t- sign((5M2) + cos^ 0^ 



(III.63) 



where we have used eqn. (IIIL53|I for the case when the momentum dependent contribution is much larger than the 
asymmetries. The in-medium correction to the group velocity being proportional to 



M^ 



10 



-21 f T 



VMeV, 
is negligible in the temperature regime in which the calculation is reliable, namely for T <^ M^. 



(III.64) 
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IV. TIME SCALES IN THE RESONANCE REGIME 

There are several important time scales that impact on the dynamics of wave-packets in the medium as revealed by 
the discussions above, but also there are two more relevant time scales that are pertinent to a plasma in an expanding 
Universe: the Hubble time scale tn ~ ^/H which is the cooling time scale T{t)/T{t) and the coUisional relaxation 
time scale trei = 1/r with T the weak interaction collision rate. Neither tn nor trei has been input explicitly in the 
calculations above which assumed a medium in equilibrium and considered self-energy corrections only up to 0{Gp). 
The damping factor that leads to the decoherence from neutral and charged current interactions has been studied in 
detail in references \Wi 122. |24L l25l| and we take this input from these references in order to compare this time scale for 
damping and decoherence to the time scales for the space-time evolution of the wave packets obtained above at the 
one-loop level. 

In the temperature regime 1 MeV <T< 100 MeV the Hubble time scale islBOJ 

t,^0.6(^) \ (IV.l) 



and the coUisional rate is estimated to bep!ol 122. 124 125| 






r ~ 0.25 G^r^~ 0.25 X 10-"" (^^^1 MeV^i,ei~1.6(^^j-^l s (IV.2) 

In order to determine the relevant time scales an estimate of the momentum spread of the initial wavepacket a is 
needed. For example, for neutrinos in the LSND experiment, the momentum spread of the stopped muon is estimated 
to be about O.OlMeV 35]. An estimate of the momentum spread in the medium can be the inverse of the mean-free 
path of the charged lepton associated with the neutrino ^|. This mean free path is determined by the electromagnetic 
interaction, in particular large angle scattering, which can be simply estimated from one-photon exchange to be 
Am/ ~ (ctg^T)-^. This estimate yields 

a ^ al^T ^ 10-^ (^^^ (MeV) . (IV.3) 

For neutrinos in the neutrinosphere of a core-collapse supernovae, the estimate for a is also ^ 10"^ MeV[3|. We will 
take a value a ^ 10^^ MeV in the middle of this range as representative to obtain order of magnitude estimates for 
the time scales, but it is straightforward to modify the estimates if alternative values of a can be reliably established. 

We now consider the large mixing angle (LMA) case to provide an estimate of the different time scales, but 
a similar analysis holds for the case of small vacuum mixing (SMA) by an appropriate change of ko^T. Taking 
fco ~ T - 3.6 MeV, a ~ lO'^MeV, \5M'^\ ~ 8 x 10-'^(eV)2, M ~ 0.25el^ we obtain the following time scales near the 
resonance region: 

i) oscillation time scales: 



47rfco 



3.8x10-^5 ; Trned^^^ (IV.4) 



\oM'^\ sn\20 

ii) dispersion and coherence time scales: 

T_L - — ^ -1.2x IQ-^^s ; Til ^- 0.25 s (IV.5) 

2^^' 2^2 Af' 

^-^-^-V^\^^-'^^' ' *"» = ^^ (IV.6) 

(T|dM^| sm26' 

iii) expansion and coUisional relaxation time scales: 

tff- 4.6x10-^5 ; i^ei -- 2.8 X 10"^ s (IV.7) 

For small vacuum mixing angle (^13) the above results are modified by taking fco ~ T ~ 7 MeV. 
In the resonance region the in-medium coherence time scale is of the same order as the Hubble time (for LMA) or 
much longer (for SMA) and and there is a large temperature variation during the coherence time scale. However, the 
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decoherence of the wave packets occurs on much shorter time scales determined by the colhsional relaxation scale and 
the coherence time scale is not the relevant one in the medium near the resonance. 

Decreasing the momentum spread of the initial wave packet a increases the dispersion and coherence time scales, 
with the dispersion scales increasing faster. The medium effects are manifest in an increase in the oscillation and 
the coherence time scales by a factor 1/ sin 20. This effect is more pronounced for 1 — 3 mixing because of a much 
smaller mixing angle. It is clear from the comparison between the coherence time scale in the medium tc,m and the 
relaxational (colhsional) time scale trei that unless a is substantially larger than the estimate above, by at least one 
order of magnitude in the case of 1 — 2 mixing, or even more for 1 — 3 mixing, collisions via neutral and charge currents 
is the main source of decoherence between the separating wave-packets near the resonance. However, increasing a will 
decrease the transverse dispersion time scale thus leading to greater suppression of the amplitude of the wave packets 
through dispersion. Furthermore, for large mixing angle sin 26 ^ 1 the oscillation scale is shorter than the colhsional 
decoherence time scale via the weak interactions trei, therefore allowing several oscillations before the wave packets 
decohere, and because the oscillation scale is much smaller than the Hubble scale the evolution is adiabatic over the 
scale trei- But for small mixing angle the opposite situation results and the transition probability is suppressed by 
colhsional decoherence, furthermore for small enough mixing angle there is a breakdown of adiabaticity. However, 
the strongest suppression of the survival Ve^e as well as the transition Ve^^ probabilities (equally) is the transverse 
dispersion of the wave packets, on a time scale many orders of magnitude shorter than the decoherence tc.m and 
the colhsional trei time scales. Unless a^ is within the same order of magnitude of |(5M^| the transverse dispersion 
occurs on time scales much faster than any of the other relevant time scales and the amplitude of the wave-packets is 
suppressed well before any oscillations or decoherence by any other process can occur. Clearly a better understanding 
of the initial momentum spread is necessary for a full assessment of the oscillation probability in the medium. 



V. DISCUSSIONS AND CONCLUSIONS 

In this article we implemented a non-equilibrium quantum field theory method that allows to study the space- 
time propagation of neutrino wave-packets directly from the effective Dirac equation in the medium. The space- 
time evolution is studied as an initial value problem with the full density matrix via linear response. The method 
systematically allows to obtain the space-time evolution of left and right handed neutrino wave packets. 

A "flavor neutrino" wave packet evolves in time as a linear superposition of wave-packets of "exact" (quasi) particle 
states in the medium, described by the poles of the Dirac propagator in the medium. These states propagate in the 
medium with different group velocities and the slow separation between these packets causes their overlap to diminish 
leading to a loss of spatial and temporal coherence. However, the time evolution of the packets also features dispersion 
as a result of the momentum spread of the wave packets jsj. 

The space time dynamics feature a rich hierarchy of time scales that depend on the initial momentum spread of 
the wave packet: the transverse and longitudinal dispersion time scales tj^ ^ T|| which are widely separated by the 
enormous Lorentz time dilation factor « (fc/M )^ with M the average neutrino mass, and a coherence time scale tc^m 
that determines when the overlap of the wave packets becomes negligible. The dynamics also displays the phenomenon 
of "freezing of coherence" which results from the competition between the separation and spreading of the wave packet 
along the direction of motion (longitudinal). For time scales larger than tii the overlap of the wave-packets freezes, 
with a large overlap in the case when tc,m S> t\\ , which occurs for a wide range of parameters. 

We have focused on studying the space-time propagation in the temperature and energy regime in which there is 
a resonance in the mixing angle in the medium, prior to BBN[22„ 24, 25, 26]. Our main results are summarized as 
follows: 

• Both the coherence and oscillation time scales are enhanced in the medium with respect to the vacuum case by 
a factor 1/ sin 20 near the resonance, where 9 is the vacuum mixing angle. 

• There are small corrections to the oscillation formula from the wave-packet treatment, but these are suppressed 
by two powers of the ratio of the momentum spread of the initial packet to the main momentum. 

• There are also small corrections to the space-time evolution from the energy dependence of the mixing angle, 
but these are negligible near the resonance region. 

• The spreading of the wave-packet leads to the phenomenon of "freezing of coherence" which results from the 
competition between the longitudinal dispersion and coherence time scales. This phenomenon is a result of the 
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longitudinal spreading of the wave-packets "catching up" with their separation. Substantial coherence remains 
frozen for ic,m ^ ■'"n • 

• We have compared the wide range of time scales present in the early Universe when the resonance is available 
for T ^ 3.6MeV [l0l 122. l24L |25L l2a| for large mixing angle. Assuming that the initial momentum spread of the 
wave-packet is determined by the large angle scattering mean free path of charged leptons in the medium Jl|, we 
find the following hierarchy between the transverse dispersion t±, oscillation Tmed, collisional relaxation trei, 
Hubble tn, in- medium coherence tc.m and longitudinal dispersion r|| time scales respectively: for large vacuum 
mixing angle sin 20 ^ 1: 

r_i < T„ied < Uel <tH <■ tc,m < T|| (V.l) 

and for small mixing angle sin 29 ^ 1 

T±_ < trel < Trned < tff < tc,m < T|| . (V.2) 

The rapid transverse dispersion is responsible for the main suppression of both the persistent and transition 
probabilities making the amplitudes extremely small on scales much shorter than any of the other scales. Only 
a momentum spread a ^ ^/\5M^ will make the transverse dispersion time scale comparable with the oscillation 
and relaxation ones. Clearly a better assessment of the momentum spread of wave-packets in the medium is 
required to provide a more reliable estimate of the wave-packet and oscillation dynamics. 
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